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Quantum discord is a more general measure of quantum correlations than entangle- 
ment and has been proposed as a resource in certain quantum information processing 
tasks. The computation of discord is mostly confined to two-qubit systems for which an 
analytical calculational scheme is available. The utilization of quantum correlations in 
quantum information-based applications is limited by the problem of decoherence, i.e., 
the loss of coherence due to the inevitable interaction of a quantum system with its 
environment. The dynamics of quantum correlations due to decoherence may be studied 
in the Kraus operator formalism for different types of quantum channels representing 
system-environment interactions. In this review, we describe the salient features of the 
dynamics of classical and quantum correlations in a two-qubit system under Markovian 
(mcmoryless) time evolution. The two-qubit state considered is described by the reduced 
density matrix obtained from the ground state of a spin model. The models considered 
include the transverse-field XY model in one dimension, a special case of which is the 
transverse-field Ising model, and the XXZ spin chain. The quantum channels studied in- 
clude the amplitude damping, bit-flip, bit-phase-flip and phase-flip channels. The Kraus 
operator formalism is briefly introduced and the origins of different types of dynam- 
ics discussed. One can identify appropriate quantities associated with the dynamics of 
quantum correlations which provide signatures of quantum phase transitions in the spin 
models. Experimental observations of the different types of dynamics are also mentioned. 
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1. Introduction 

Interacting quantum systems are characterized by the presence of correlations be- 
tween the different parts. The correlations are of two types: classical and quantum. 
Since the early days of quantum mechanics, the idea of entanglement has been in- 
voked to probe the origins of quantum correlations^. In later years, various quanti- 
tative measures of entanglement have been proposed and the rol e of e ntanglement 
in a number of quantum information- based protocols highlighted EEE] ^ quantum 
state, pure or mixed, is either entangled or separable, the latter condition implying 
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the absence of quantum correlations as measured by entanglement in the state under 
consideration. In recent years, a new measure of quantum correlations, the quan- 
tum discord (QD), has been proposed based on the information theoretic concept 
of mutual information IHEE The QD is quantified by the difference between two 
quantum extensions of the classical mutual information. The two representations 
are identical in the classical domain. 

In classical information theory, the total correlation between two random vari- 
ables A and B is measured by their mutual informational 

I (A, B) = H(A) + H(B) - H(A, B) (1) 

The random variables A and B take on the values 'a' and respectively with prob- 
abilities given by the sets {p a } and {pb}- The probability distribution {p a b} defines 
the outcome when joint measurements are carried out. The variables A and B are 
correlated when {p a b} does not have a product form {p a x pb}. In Eq. (TT]), H{A) = 
-HaPa l°g'2Pa, H (B) = - £ 6 Pb log 2 and H(A, B) = - J2 a ,b Pab log 2 Pab are the 
Shannon entropies for the variables A, B and the joint system AB respectively. The 
probabilities p a , Pb and p a b satisfy the relations p a = J2bPab and pb = J2 a p a b- The 
Shannon entropy of a random variable A quantifies our ignorance about A before 
we measure its value or equivalently it provides a measure of how much informa- 
tion is gained on an average after a measurement is carried out An alternative 
representation of the classical mutual information is given by I^ES 

J(A,B) = H(A)-H(A\B) (2) 

where H(A\B) is the conditional entropy and quantifies our lack of knowledge of 
the value of A when that of B is known. The exact equivalence of the expressions 
in Eqs. (QJ and © can be demonstrated using the Bayes rule p a b = P a \bPb and the 
definition H(A\B) — — J2 a bP ab ^°&2Pa\b of the conditional entropy. 

The generalization of the classical mutual information to the quantum case is 
achieved by replacing the classical probability distribution and the Shannon en- 



tropy by the density matrix p and the von Neumann entropy S(p) = — tr (p\og 2 p), 
respectively. The quantum generalizations of Eqs. ([TJ and ([2]) are given by 

I {pab) = S(p A ) + S(p B ) - S(pab) (3) 
J (Pab) = S (p A ) - S (p A \p B ) (4) 



where S (pab) is the quantum joint entropy and S (pa\pb) the quantum conditional 
entropy. The latter quantity is, however, ambiguously defined as the magnitude of 
the quantum conditional entropy (ignorance of A once B is known) depends ex- 
plicitly on the type of measurement carried out on B. Since different measurement 
choices yield different results, Eqs. ((3j) and Q arc no longer identical. We consider 
von Neumann-type measurements on B defined in terms of a complete set of orthog- 
onal projectors {nf } corresponding to the set of possible outcomes i. The state of 
the system after the measurement is given by 

ft = (/®nf)p A fl(j®nf)/p j (5) 
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with 

Pi =tr((l®Ilf)p AB (l®Il?)) (6) 

/ denotes the identity operator for the subsystem A and pi is the probability of 
obtaining the outcome i . From Eq. fl|) , an alternative expression of quantum mutual 
information is given by 

j(p AB ,{llf})=S(p A )-S(p AB \{llf}) (7) 

The quantum analogue of the conditional entropy is 

s(p AB \{n*}) = J2piS(pi) (8) 

i 

Henderson and Vedral ^ have shown that the maximum of J [pab, {nf }) w.r.t. 
{nf } provides a measure of the classical correlations (CC), C [pab), i-e., 

C(p AB )= max {j{p AB ,{nf})) (9) 
{nf} 

The difference between the total correlations / {pab) (Eq. ©) and the CC, C (pab), 
defines the QD, Q (pab)- 

Q (pab) = I (pab) - C (pab) (10) 

The QD is defined for bipartite systems only and due to the computational difficulty 
of carrying out the extremization process in Eq. ©, the calculation of the QD is 
mostly confined to two-qubit systems. For any pure state, the QD reduces to the 
entropy of entanglement and the total correlations, measured by the mutual 
information, are equally divided between the classical and quantum correlations, 
i.e., C (pab) = Q (pab) = \l (pab)- 

In the case of mixed states, however, the QD and the entanglement provide dif- 
ferent measures of quantum correlations. In fact, there are mixed states which are 
separable, i.e., unentanglcd but for which QD is non-zero. There is thus a shift in 
focus in the case of the QD from the separability versus entanglement criteria to 
the issue of classical versus quantum correlations. The relationship between the QD, 
entanglement and classical correlations is not as yet clearly understood even for the 
simple two-qubit system. The QD is not alw ays larger than entanglement and is 
not always less than classical correlations ^21. The QD defined by Eqs. ([3]), ((9|) and 
(|10j) does not provide a measure of quantum correlations in multipartite systems. 
Recently, the concept of relative entropy has been utilized to obtain measures of 
classical and non-classical correlations in a given quantum state I 13 | 14 | The relative 
entropy between two quantum states x and y is given by S^H?/) = tr ^x\og 2 . It 
is a non-negative quantity and hence serves as a "distance" measure of the state x 
from the state y. The relative entropy of entanglement is the distance of the entan- 
gled state p to the closest separable state a. Similarly, the discord is measured by the 
distance between the quantum state and its closest classical state. Analogous defini- 
tions exist for quantum dissonance (non-classical correlations for separable states), 
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total mutual information and classical correlations. The measures are valid in the 

multipartite case and for arbitrary dimensions of the system. Since the different 

types of correlations have a common measure, additivity relations connecting the 

correlations can be derived . O ther alternative measu res o f non-classical correlations 

include the geometric QD an d th e Gaussian QD 021. In this review, we confine 

ourselves to the original definition of the QD as embodied in Eqs. ([3]), ([9]) and (fT0|) . 

Condensed matter systems like molecular magnets (represented by spin clusters) 

and spin chains have been extensively studied to characterize as well as quantify the 

quantum correlations present in the ground and thermal states of the spin systems 

| 2 | 3 | 17 | 18 ] while the bulk of the studies are devoted to entanglement in its various 

for ms, there are now several studies on the QD pro perties of well-known spin mod- 
els | 7|17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27|28 | 29 | 30 | 3ir32 | Many of thcse studies focus on 

how the QD and associated quantities provide signatures of quantum phase tran- 
sitions in specific spin systems. Quantum phase transitions (QPTs) in interacting 
systems occur at T = and are brought about by tuning a non-thermal parame- 
ter g, e. g., pre ssure, chemical composition or external magnetic field, to a special 
value g e I 33 | 34 | -p^g ^ rans j^i on j s driven by quantum fluctuations and brings about 
qualitative changes in the ground state wave function at the transition point. It 
is then reasonable to expect that the quantum correlations present in the ground 
state would provide signatures of the occurrence of a QPT. Apart from the effect of 
changing parameters on the quantum correlations of a system, the inevitable inter- 
action of the system with its environment results in decoherence, i.e., a destruction 
of quantum properties including correlations The dynamics of entanglement and 
QD under >ystoin environment interaction,- have been investigated in a number of 
recent studies | 20 | 29 | 35 | 36 | 37 | 38 | 39 | 40 | 41 | Q ne f eature w hi c h emerges out of such 

studies is that the QD is more robust than entanglement in the case of Markovian 
(mcmoryless) time evolution. The dynamics may bring about the sudden disap- 
pearance of entanglement at a finite time termed the 'entanglement sudden death' 
| 35 | 36 | The QD, on the other hand, decays in time but vanishes only asymptot- 
ically | 20 | 29 | 37 | 39 | 40 | 41 ] Also, under Markovian time evolution and for a class of 
initial states , the d ecay rates of the classical and quantum correlations exhibit sud- 
den changes Mo | 39 | In two recent studies, the dynamics of the mutual information 
I (pab), the classical correlations C (pab) and the quantum correlations Q (pab), 
as measured by the QD, have been studied in two-qubit states the density matri- 
ces of which are the reduced density matrices obtained from the ground states of 
the transverse- field Is ing model (TIM) and the transverse-field XY model in one 
dimension (Id) I 29 | 40 | time evolution brought about by system-environment 
interactions is assumed to be Markovian in nature and the quantum channels, rep- 
resenting qubit-environment interactions, include amplitude damping (AD), bit- flip 
(BF), phase- flip (PF) and bit-phase-flip (BPF). A significant outcome of the studies 
is the identification of appropriate quantities associated with the dynamics of the 
correlations which signal the occurrence of a QPT. The TIM in Id is a special case 
of the transverse-field XY model in Id. The latter model has a rich phase diagram 
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exhibiting QPTs | 42 | 43 | 44 | 45 | 46 | 47 ] j n f ac ^ k th the spin models are well-known 
statistical mechanical models which illustrate QPTs. In this review, we combine the 
themes of quantum correlations, decoherence and QPTs in the study of spin models 
like the TIM in Id, the transverse-field XY chain and the XX Z spin chain. In Sec- 
tion 2, we describe the spin models and the QPTs associated with them. In Section 
3, we introduce the Kraus operator formalism for describing the time evolution of 
open quantum systems, i.e., systems interacting with specific environments. Some 
representative quantum channels are described as also the dynamics of mutual in- 
formation, clas sical a nd quantum correlations. Section 4 describes the main results 
obtained so far I29 | 40 ] as we rj ag some new re sults. The methodology reported in the 
review is general in nature and applicable to other condensed matter models ex- 
hibiting QPTs. In Section 5, some concluding remarks are made and future research 
directions pointed out. 

2. Quantum Phase Transitions in Spin Models 
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Fig. 1. Criticality of the transverse-field XY model in the (7, A) parameter space. The transverse- 
field Ising model (7 = ±1) has a QCP at A c = 1. The line A = 1 represents criticality of the XY 
model with universality class the same as that of the transverse-field Ising model. The 7 = line 
(XX model) represents the anisotropy transition line for A S [0, 1] with the transition belonging 
to a new universality class. 



The fully anisotropic Hcisenberg spin chain in a magnetic field is described by 
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the Hamiltonian 

L L 
i=l i=l 

where of {a = x, y,z) is the Pauli spin operator at the ith site, J x , J y and J z 
are the strengths of the nearest-neighbour (n.n.) exchange interactions and h the 
external magnetic field. L denotes the total number of spins in the chain. A num- 
ber of spin models are special cases of the fully anisotropic model: (i) XX Z spin 
model in a magnetic field (J x = J y ^= J z , h ^ 0), (ii) transverse-field XY model 
(J x 7^ J y , J z = 0, h ^ 0), (Hi) transverse-field XX model (J x = J y , J z = 0, h ^ 0) 
and (iv) the TIM (J y = J z = 0, h ^ 0). One can also consider cases for which h = 0. 
The TIM in this case reduces to the Ising model which has no quantum character. 

The transvcrsc-ficld XY model in Id describes an interacting spin system for 
which many exact resu lts on the ground and excited state properties including spin 
correlations are known B2 | 43 | 44 | correS p nding Hamiltonian is written as 

H XY = ~ / Z{( 1 + 7)«i + (1 - 7)*M + i} - E < ( 12 ) 

i=l i=l 

where 7 is the degree of anisotropy (— 1 < 7 < 1) and A is inversely proportional 
to the strength of the transverse magnetic field in the z direction (A > 0). The 
Hamiltonian satisfies periodic boundary condition, i.e., L + 1 = 1 and is transla- 
tionally invariant. Two special cases of the XY model are the TIM with 7 = ±1 
and the isotropic XX model (7 = 0) in a transverse magnetic field. For the full 
range of values of the anisotropi c parame ter, Hxy can be diagonalizcd exactly in 
the thermodynamic limit L — > 00 1 42 | 43 | 44 [ ^^jg j s achieved via the successive appli- 
cations of the Jordan- Wigner and Bogoliubov transformations. Figure [T] describes 
the critical transitions associated with the transverse- field XY model in the (7, A) 
parameter space. For non-zero values of 7, a second-order QPT occurs at the critical 
point A c = 1 separating a ferromagnetic ordered phase (A > 1) from a quantum 
paramagnetic phase (A < 1). The transition is characterized by the order param- 
eter {a x ), the magnetization in the x direction, which has a non-zero expectation 
value only in the ordered ferromagnetic phase (A > 1). The magnetization in the z 
direction, (er 2 ), is non-zero for all values of A with its first derivative exhibiting a 
singularity at the critical point A c = 1. When < I7I < 1, the critical point transi- 
tion occurring at A c = 1 belongs to the Ising universality class. For A € [0, 1], there 
is another QPT, termed the anisotropy transition, at the critical point 7 = 0. The 
transition belongs to a different universality class and sep arates two ferr omagnetic 
phases with orderings in the x and y directions respectively | 42 | 43 | 44 | 45 | 45 ] pj]y[ 
Hamiltonian in Id is obtained from Eq. (|12[) by putting 7 = 1 and is given by 



L 

iW = -A5>f< +1 -5>? (13) 
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When the parameter A = 0, all the spins are oriented in the positive z direction in 
the ground state whereas in the extreme limit A = oo the ground state is doubly 
degenerate with all the spins pointing in either the positive or the negative x direc- 
tion. As mentioned before, a QPT occurs at the critical point A c = 1 with (a x ) ^ 
in the ordered ferromagnetic phase (A > 1). The XXZ spin chain in a zero magnetic 
field has the Hamiltonian 



Two QPTs, one first-order and the other infinite-order, occur at the respective 
points A = 1 (ferromagnetic point) and A = — 1 (antifcrromagnctic point). At 
the latter point, the XXZ spin chain undergoes a QPT from an antifcrromagnctic 
phase (A < —1) to an XY phase for — 1 < A < 1. 

In recent yea rs, quantum informatio n-related measures like entanglement, dis- 
cord and fidelity l ^ | 3 | 9 | 19 | 48 | 49 | 50 | 51 | 52 | ^ been shown to provide signatures of 
QPTs. An nth order QPT is characterized by a discontinuity/divergence in the 
nth derivative of the ground state energy with respect to the tuning parameter g 
as g — > g c , defining the transition point. At a first-order QPT point, a ppropriate 
entanglement measures and the QD are known to become discontinuous I ' ' ' ' 
whereas a second-order (critical point) QPT is signaled by the discontinuity or di- 
vergence of the first derivative of the quant um corre lation measures with respect 
to the tuning parameter at the critical point I 9 | 19 | 22 [ a q uan t U m critical point, 
quantum fluctuations occur on all length scales leading to a divergent correlation 
length. The ground state energy and related quantities become non-analytic as the 
tuning parameter g tends to the critical point g c . The influence of a QPT extends 
into the finite temperature part of the phase diagram so that experimental detection 
of the QPT is possible. 

We now quote some major results on the signatures of QPTs provided by quan- 
tum correlation measures like entanglement and the QD in the cases of spin models 
like the TIM in Id, transverse- field XY and XXZ spin chains. Fo r more detailed in- 
formation, the reader is referred to the original papers and reviews ^^^^^'^ . In 
the case of the TIM in Id, the n.n. concurrence, a measure of pairwisc entanglement, 
reaches its maximum value close to the QCP A c = 1. The ncxt-ncarcst-ncighbour 
(n.n.n.) concurrence, on the other hand, has its maximum value at A c = 1- The 
first derivative of the n.n. concurrence w.r.t. the tuning parameter A exhibits a 
logarithmic divergence as A — > A c . The pairwisc entanglement as measured by the 
concurrence does not become long-ranged as A — > A c , in fact it falls to zero beyond 
the n.n.n. distance. Both the n.n. and the n.n.n. QD attain their maximal values 
close to the critical point A c . The first derivative of the n.n. (n.n.n.) QD with respect 
to A becomes discontinuous (singular) at the critical point A c = 1. The magnitude of 
the classical correlations has a monotonic dependence on A starting with low values 
for small A. The transverse-field XY chain belongs to the universality class of the 
TIM in Id away from the isotropic limit 7 = 0. In the case of the zero-field XXZ 




(14) 



; = 1 
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chain, the n.n. concurrence is maximal at the infinite-order critical point A c = — 1. 
The n.n. QD is found to be maximal (with a discontinuity) and the classical cor- 
relations become minimal (with a kink) at A c = — 1. At the first-order transi tion 
point A = 1, the n.n. classical and quantum correlations arc discontinuous SHU p or 
the transverse-field XY chain, Maziero et. al. derived analytic expressions for 
classical and quantum correlations for spin-pairs separated by an arbitrary distance 
for both temperatures T = and T^O. The T = QD for spin pairs beyond the 
n. n.n. distance is able to signal a QPT whereas pairwise entanglement for the same 
spin pairs is unable to do so. Pairwise entanglement is typically short-ranged even 
close to criticality. Spin chains like the transverse-field XY chain and the XXZ 
chain in the presence of domain walls are characterized by a long-range decay of the 
QD as a function of the spin-spin distance close to the quantum critical point ^31. 
The decay rate of the QD has a noticeable change as the critical point is crossed. In 
a recent study ^3, the thermal QD is shown to be maximal and its first derivative 
with respect to the tuning parameter discontinuous at the quantum critical point 
A = -1 of the XXZ chain for both T = and T > 0. The entanglement of for- 
mation is maximal at the critical point only for T = 0, the maximum shifts from 
A = -1 for T > 0. The thermal QD also signals the QPT at A = 1. Thus, the QD 
has the important property of being able to detect the critical points of QPTs at 
finite temperatures. In contrast, both entanglement and thermodynamic quantities 
fail to signal QPTs for T > 0. The special property of the QD could be useful in 
the experimental detection of critical points. 



3. Dynamics of Correlations 

We next consider the interaction of the chain of qubits (spins) with an environment. 
The initial state of the whole system at time t = is assumed to be of the product 
form, i.e., 

p(0) = p s (0)® Pe (0) (15) 

with the density matrices p s and p e corresponding to the system (spin chain) and 
the environment respectively. We assume that the environment is described in terms 
of L independent reservoirs each of which interacts locally with a qubit constituting 
the spin chain. The two-qubit reduced density matrices, p rs and p re , are obtained by 
taking partial traces on p s and p e respectively over the states of all the qubits other 
than the two chosen qubits. The two-qubit reduced density matrix p r (0) obtained 
from Eq. ()15() can be written as 

Pr (0) = p rs (0)®p re (0) (16) 

The quantum channel describing the interaction between a qubit and its loca l envi- 
ronment can be of various types: AD, phase damping, BF, PF, BPF etcJ ^M With 
the initial reduced density matrix given in Eq. (fTB"|) . the objective is to determine 
the dynamics of the two-qubit classical and quantum correlations (in the form of 
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the QD) under the influence of various quantum channels. The time evolution of 
the closed quantum system consisting of both the system and the environment is 
given by 

Pse (t) = U(t) Pse (0)U\t) (17) 

where U(t) is the unitary evolution operator generated by the total Hamiltonian H 
(U = e ~ tHt / h ^ . H is written as H = H s + H e + H se where H s and H e are the bare 
system and the environment Hamiltonians respectively and H se the Hamiltonian 
describing the interactions between the system and the environment. The time 
evolution of the system s under the influence of the environment e is obtained by 
carrying out a partial trace on p se (t) (Eq. (|17[1 ) over the environment states, i.e., 

p s (t) = Tr e [U(t)p se (0)U\t)] (18) 

In Eq. (JTHJ) , p S e{0) = p(0) from Eq. (fT5| . Let \ek) be an orthogonal basis spanning 
the finite-dimensional state space of the environment. With the initial state of the 
whole system given by Eq. (fT5)l , 

p s [t) = 5^< ejb |l7 0,(0) ® p e (0)\ tf\e k ) (19) 
fe 




Fig. 2. BPF and BF channel: Decay of mutual information J (pab) (solid line), classical cor- 
relations C (pab) (dashed line), and the QD Q (pab) (dot-dashed line) as a function of the 
parametrized time p = 1 — e~ et for A = 0.7,7 = 0-7 (BPF channel) and A = 0.7,7 = —0.7 
(BF channel). Also, p BC = 0.114. (inset) Decay of mutual information / (pab) (solid line), classical 
correlations C {pab) (dashed line) and quantum discord Q {pab) (dot-dashed line) as a function 
of the parametrized time p = 1 - e~ Bt for A = 0.7, 7 = -0.7 (BPF channel) and A = 0.7, 7 = 0.7 
(BF channel). 

Let p e (0) = |eo)(eo| be the initial state of the environment. Then 

p B {t) = ^E k p a (p)El (20) 

k 
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where Ek = (ek\U \eo) is the Kraus operator which acts on the state space of the 
system only 121251 Let i = 1,2, <i, define the basis in the state space of the 

syste m s. There are then at most d 2 independent Kraus operators Ek-, k = 0, d 2 — 1 



systo 
8 53 



The unitary evolution of s + e is given by the map: 

|0i)|e o ) -> E o |0i)|e o ) + ...+E d 2_ 1 \ ( j) 1 )\e d 2_ 1 ) 
|02>|e o ) -> ^o|02>|e o > + + E d 2_ 1 \ ( j )2 )\e d 2_ 1 ) 



\<t> d )\eo) ->■ ^o|0d)|eo) + ... + ^-i|0d)|e d 2„ 1 ) (21) 
In compact notation, the map is given by 

U\4H)\eo) = ^2E k \4H)\e k ), i = l,2,...,d (22) 

In the case of N system parts with each part interacting with a local independent 
environment, Eq. (|20p becomes 

(23) 

where E^ is the /c a th Kraus operator with the environment acting on the system 
part a. The specific form for p s (t) arises as the total evolution operator can be 
written as U(t) = U\{t) <£> J7 2 (t) <g> ... <g> J7jv(t). Following the general formalism of 
the Kraus operator repre sentat ion, an initial state, p rs (0), of the two-qubit reduced 
density matrix evolves as^2~ZI 

p rs {t) = Y J E li . v Prs(Q)El,u (24) 

where the Kraus operators E^ = E^ ® E v satisfy the completeness relation 
y t v E^EJt v = I for all t. We now briefly describe the various quantum chan- 
nels and write down the corresponding Kraus operators. A fuller description can be 
obtained from Refs. [8l 137). 

(i) AD Channel. The channel describes the dissipative interaction between a system 
and its environment resulting in an exchange of energy between s and e so that s is 
ultimately in thermal equilibrium with e. The s + e time evolution is given by the 
unitary transformation 

|0>.|0) e |0) s |0) e (25) 
|l) s |0) e -> v^|l) s |0) e + v/p|0),|l) e (26) 

where |0) s and |l) s are the ground and excited qubit states and |0) e , |l) e denote 
states of the environment with no excitation (vacuum state) and one excitation 
respectively. Eq. (f25"T) stipulates that there is no dynamic evolution if the system 
and the environment are in their ground states. Eq. (|26|) states that if the system 
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qubit is in the excited state, the probability to remain in the same state is q and 
the probability of decaying to the ground state is p (p + q = 1). The decay of 
the qubit state is accompanied by a transition of the environment to a state with 
one excitation. The qubit states may be two atomic states with the excited state 
decaying to the ground state by emitting a photon. The environment on acquiring 
the photon is no longer in the vacuum state. With a knowledge of the map equations 
(Eqs. (|2"5|) and (j2"?5)) ). the Kraus operators for the AD channel can be written as 

where q = 1 — p. The Kraus operators for the two distinct environments (one for 
each qubit) have identical forms. In the case of Markovian time evolution, p is given 
by p = 1 — e~ et with 9 denoting the decay rate. 

(ii) Phase Damping (dephasing) Channel. The channel describes the loss of quantum 
coherence without loss of energy. The Kraus operators are: 



with q = 1 — p and p = 1 



-01 




Fig. 3. Variation of p sc as a function of A with 7 = 0.5 for the BPF channel, (inset) The first 
derivative of p sc w.r.t. A diverges as the QCP A c = 1 is approached. 



(iii) BF, PF and BPF channels. The channels destroy the information contained in 
the phase relations without involving an exchange of energy. The Kraus operators 
are 

E = ^ (J J)l E 1 = ^p~j2a i (29) 

where i = x for the BF, i = y for the BPF and i = z for the PF channel with 
q' = 1 — p/2 and p = 1 — e~ et . The expanded forms of the Kraus operators are: 
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BF: 



PF: 



BPF: 



E Q = 
Ei = 

E = 
Ei = 

E = 
Ei = 



o yi^/s 
o vWs 

VW2 o 

o 

VW2 o 

-Vp/2 



(30) 



(31) 







^T~pj2 

-i^/p/2 
is/pf2 



(32) 



As shown in Ref. [8], the phase damping quantum operation is identical to that of 
the PF channel so that we will consider only one of these, the PF channel, in the 
following. 

In the case of the transverse-field XY chain, including the TIM, the reduced 
density matrix pab has the form of X-states with 

/a00/\ 
ObzO 

P ab= QzbQ (33) 
\f00dj 

For the XX Z spin chain, the element / = 0. In Eq. (|33[) . A, B represent the two 
individual qubits and z, / are real numbers. The eigenvalues of pab are^ 



with 



A = \ |(1 + 03) + ^cl + (ci-c 2 f 
Ai = \ |(1 + c 3 ) - ^4cl + (ci-c 2 y 
A 2 = j (1 - c 3 + ci + c 2 ) 
A3 = - (1 - c 3 - c x - c 2 ) 



ci = 2z + 2/ 
c 2 = 2z - 2/ 
C3 = a + d — 2b 
C4 = a — d 



(34) 



(35) 
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The mutual information / (pab) can be written as^ 

3 

/ (pab) = S( PA ) + S{p B ) + Y, X » 1o 8'2 A Q (36) 

where 

-5 (pa) = S (p B ) = — log 2 — — log 2 — ^— (37) 

With expressions for / (pab) and C (pab) given in Eqs. (j36]) . (|37)l and ©, the QD, 
Q (pab), (Eq. ((10)0 can in principle be computed. The difficulty lies in carrying out 
the maximization procedure needed for the computation of C (pab)- It is possible 
to do so analytically when pab is o f the form given in Eq. (|33|) resulting in the 
following expressions for the QD 



Q( P ab) =min{Q 1 ,Q 2 } (38) 



where 



S (p B ) - S (pab) - alog 2 ~r~r - folo 8'2 — ~r - dl °g2 XTT ~ &lo &2 



a + b a + b d + b d + b 

(39) 

and 

Q 2 = S ( PB ) - S (pab ) - A+ log 2 A+ - A_ log 2 A_ (40) 

with A± = | (1 ± r) and r 2 = (a - df + 4 (\z\ + |/|) 2 

Let pij be the reduced density matrix for two spins located at the site i and 
j respectively, p^ has the form shown in Eq. ()33|) . The matrix elements of p^ 
can be expressed in terms of single-site magnetization and two-site spin correlation 
functions. In the case of the transverse-field XY spin chain, the matrix elements 



d 



4 2 4 

1 (O 
4 2 4 

* = \ {{of ^) + 

/ = J «of *?) - (af a})) (41) 

The sing le-site magnetization (er 2 ) in the case of the transverse-field XY model is 
given by I 46 | 19 | 



7T 



(42) 
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where 



= J (7Asin(/)) 2 + (1 + Acos< 



(43) 



describes the energy spectrum. The spin- spin corre lation functions are obtained 
from the determinants of Toeplitz matrices | 19 | 43 | 44 | ag 



where 



G — i G_2 • ' • G— r 

Go G—\ ■ ■ ■ G_ r +i 



Gy Gy—\ • ■ ■ G\ 
(ofaf +r ) = (<r z ) 2 - G r G- r 



C_J J - 


-2 G r _3 


■■■ GLi 




Go • • 


G- r +2 


G 2 


G 1 ■■ 


G-r+3 



n 1 C aa i ^ (1 + Acos^) 7 A r 
G r = — d(pcos{r(j)) 

71" Jo w 



dtp sa\{r<f) 

7T ,/ W 



(44) 



(45) 



with r = \i — j\ being the distance between the two spins at the sites i and j (for 
the nearest- neighbour case, r = 1). 



0.18 




0.2 0.4 0.6 0.8 

P 



Fig. 4. PF channel: Decay of mutual information I (pab) (solid line), classical correlations 
C (pab) (dashed line), and the QD Q (pab) (dot-dashed line) as a function of the parametrized 
time p = 1 — e~ 0t for A = 0.5 and 7 = 1. Q (pab) is greater than C (pab) hi the parametrized 
time interval Ap. 
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4. Results 

Three general typ es o f dynamics under the effect of decoherence have been observed 
in an earlier study ^Sl (%) C (pab) is constant as a function of time and Q (pab) de- 
cays monotonically, (ii) C {pab) decays monotonically over time till a parametrized 
time p sc is reached and remains constant thereafter. At p sc , Q (pab) has an abrupt 
change in the decay rate which has been demonstrated in actual experiments 1 5 5 1 5 6]_ 
Also, a parametrized time interval exists in which Q (pab) has a magnitude greater 
than that of C (pab) and (Hi) both C (pab) and Q (pab) decay monotonically. 
Mazzola et al. have obtained the significant result that under Markovian dy- 
namics and for a class of initial states the QD remains constant in a finite time 
interval < t < t. In this time interval, the classical correlations C (pab) decay 
monotonically. Beyond t = t,C (pab) becomes constant whereas the QD decreases 
monotonically as a function of time. In the case of the spin models under consid- 
eration, the rules of evolution for the coefficients arc obtained from Eqs. (|3l?l) . (|41[) 
and (f2"4"| with the choice of the Kraus operators dictated by the specific type of 
quantum channel (see Eqs. ([3U)l . (f3"Tj) and (f3"2")) ). The evolution rules are given by 
BPF Channel: 



BF Channel: 



PF Channel: 



d(p) = ci(0)(l-p) 2 




02 (p) = c 2 (0) 




c 3 (p) = c 3 (0)(l - p) 2 




C4(p) = C4(0)(l-p) 


(46) 


ci(p) =ci(0) 




c 2 (p) = c 2 (0)(l - p) 2 




c 3 (p) = c 3 (0)(l - p) 2 




C4(p) = C4(0)(l-p) 


(47) 


Cl (p) = Cl (0)(l - p) 2 




ca(p) = c 2 (0)(l- P ) 2 




C3(p) = c 3 (0) 




c 4 (p) = c 4 (0) 


(48) 



For the transverse-field XY model, the dynamical evolution of the mutual infor- 
mation / (pab), the classical correlations C (pab) and the QD Q (pab) in the case of 
the BPF channel are shown as a function of the parametrized time p (p = 1 — e~ et ) 
in figure [2] for A = 7 = 0.7. The dynamics are similar to type (ii) where sudden 
changes in the decay rates of the classical correlations and the QD occur at a 
parametrized time instant p sc . The inset of the figure [2] shows the time evolution 
of / (pab), C (pab) and Q (pab) for A = 0.7, 7 = —0.7 which are described by the 
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Fig. 5. Variation of p sc as a function of A with 7 = 1 for the PF channel, (inset) The first 
derivative of p sc w.r.t. A diverges as the QCP A c = 1 is approached. 



type (Hi) dynamics, i.e. both C (pab) and Q (pab) decay monotonically. Q (pab) 
tends to zero in the asymptotic limit p — > 1 whereas C (pab) and / (pab) have finite 
values in the same limit. The dynamics in the case of the BF channel have an inter- 
esting correspondence with the dynamics of the BPF channel. Type (ii) and Type 
(Hi) dynamics are obtained in the case of the BPF (BF) channel for +vc (-ve) and 
-ve (+ve) values of 7 respectively ^1. An analytical expression of the parametrized 
time p sc at which the sudden changes in the decay rate of C (pab) and Q(pab) 
take place can be obtained for the BPF/BF channels from Eqs. (|46|) and (|47|) re- 
spectively®!. The dynamics of / (pab), C (pab) and Q (pab) in the case of the PF 
channel bel ong to type (ii) with Q (pab) > C (pab) during a parametrized time 
interval SOI. The classical correlations remain constant at a value I (pab) p= i, the 
mutual information of the fully decohered state, in the parametrized time interval 
p sc < p < 1. In the interval < p < p sc , the classical correlations decay with time. 
On the other hand, the quantum correlation, as measured by the QD, undergoes a 
sudden change in its decay rate at p = p sc and goes to zero in the asymptotic limit 
p — > 1. The type of dynamics in the case of the PF channel remains the same under 
a change in sign of the anisotropy parameter 7. 

Figure [3] shows the plots of p sc and the first derivative of p sc w.r.t. A (inset) 
as functions of A for the BPF channel (7 > 0). The first derivative of p sc w.r.t. A 
diverges as the QCP A c = 1 is approached indicating a QPT. In the case of the 
PF channel, the dynamics are of type (ii) and Q (pab) is greater than C (pab) in 
a parametrized time interval. Also, the first derivative of p sc w.r.t. A for a fixed 
value of th e anisotropy parameter 7 (both +vc and — ve) diverges at the QPT point 
A c = 1 40. On the other hand, the first derivative of p sc w.r.t. 7 for a fixed value of 
A shows a discontinuity at the anisotropy tr ans ition point 7 = indicating a QPT 
for all the three channels BPF, BF and PF HS. In the case of the PF channel, p sc 
has a symmetric variation across the anisotropy transition point 7 = whereas a 
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parametric time p sc exists only for 7 > (7 < 0) in the case of the BPF (BF) 
channel. A fuller description and discussion of the results on the dynamics of the 
classical and quantum correlations are given in Ref. [40] . In the case of the AD 
channel, the dynamics are of type (Hi) with all the three correlations, I (pab), 
C (pab) and Q (pab) decaying asymptotically (p —> 1) to zero. In the case of the 
TIM in Id, the signature of quantum criticality at A c = 1 is obtained only in the 
cases of the BPF and PF channels. In the case of the PF channel, the dynamics are 
of type (ii) (figure [4]). Figure [5] shows the variation of p sc as a function of the tuning 
parameter A. The inset of the figure shows that the first derivative of p sc w.r.t. A 
exhibits a divergence as the quantum critical point A c = 1 is approached. 



0.9 
0.8 
0.7 
0.6 



0.4 
0.3 
0.2 
0.1 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

A 

Fig. 6. Variation of p sc as a function of A for the BPF/BF (A < —1) (line with solid squares) 
and the PF (— 1 < A < 1) (line with solid circles) channel, (inset) In the case of the PF channel, 
the first derivative of p BC w.r.t. A diverges as the ferromagnetic point A = 1 is approached. 

In the case of the XXZ spin chain, the critical point analysis in the presence 
of decoherence is carried out utilizing the results obtained in Refs. [HI EE] ■ Figure [S] 
exhibits the variation of p sc as a function of the anisotropy constant A (Eq. (|14[) ) 
in the cases of the BPF/BF (solid squares) and the PF (solid circles) channels. In 
the case of the PF channel, the first derivative ^j^- diverges as the ferromagnetic 
point A = 1 is approached. One also notes that p sc has a finite value when A is 
< -1 (-1 < A < 1) in the case of the BF/BPF (PF) channel. 

5. Concluding Remarks 

In this review, we have considered three spin models, namely, the TIM in Id, the 
transverse-field XY and the XXZ spin chains. The two-spin reduced density ma- 
trix obtained from the ground state in each case has the general form of Af -states 
(further simplified in the case of the XXZ chain). The matrix elements of the re- 
duced density matrix can be expressed in terms of single-spin expectation values 
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and two-spin correlation functions which are mostly known. Analytical formulae 
are available to compute the mutual information / (pab), the classical correlations 
C (par ) and the Q D Q (pab) in the two-spin state described by the reduced density 
matrix | 9 | 10 | H |12] j n t nc presence of system-environment interactions resulting in 
decohcrence, the time evolution o f the reduced density matrix can be determined 
in the Kraus operator formalism The S pi n models under consideration ex- 

hibit QPTs at specific values of a tuning parameter, e.g., magnetic field strength 
or an anisotropy parameter. While there is a large num ber of studies on how quan- 
tum correlation measures provide signatures of QPTs l ^ l ^ | 17 | 18 |50| there j s little 
inves tigatio n so far on appropriate indicators of QPTs in the presence of decohcr- 
ence l ^9 | 40 [ j n th e case f Markovian time evolution, quantities like ^jjf- and 
(transvcrsc-ficld XY model), associated with type (ii) dynamics, diverge/become 
discontinuous as a quantum critical point is approached. Dynamics similar to type 
(ii) are associated with specific quantum channels, e.g., in the case of the AD 
channel the dynamics are of type (Hi) which do not provide any indication of the 
occurrence of a QPT. There are several directions in which the investigations re- 
ported in Rcfs. [29l 00] can be extended. The analysis has so far been carried out 
for the A"-states but si nce an analytic computational scheme for a general two-qubit 
state is now available HJ one could probe the dynamics of classical and quantum 
correlations in general two-spin states and look for signatures of QPTs. In Ref. [29] . 
the reduced density matrix describes n.n. spin pairs and the analysis is carried out 
at T = 0. Extension to further-neighbours and finite temperature cases has been 
carried out in Ref. [3D] , The more general case of non-Markovian time evolution 
vis-a-vis signatures of QPTs is yet to be addressed. The dynamics of multipartite 
quantum correlations constitute another important ar ea o f study. Models of inter- 
acting fermions and bosons exhibit a variety of QPTs ^3 w hich could be signaled 
by quantities ass ociate d with the dynamics of quantum correlations. Finally, re- 
cent experiments I 55 | 56 | on dynamics of quantum correlations set new challenges 
in exploring the interconnections between quantum correlations, decoherence and 
quantum phase transitions. 
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